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A calculation of the next-to-leading order cross section for the inclusive hadropro- 
duction of charm quarks including heavy quark mass effects in the hard scattering 
cross sections is presented. It is described how the massive hard scattering cross 
sections are constructed from the corresponding cross sections in a fixed order 
calculation where coUinear pieces associated with the heavy quark are not yet 
I , subtracted. By adjusting suitable subtraction terms a massive theory with MS 

^ 1 subtraction is established which approaches the massless theory with increasing 

transverse momentum. These results will allow to have a more solid comparison 
with recent data for the inclusive D*^ cross section in pp collisions from the CDF 
' collaboration at the Tevatron at y/S = 1.96 TeV. 

oo 

o 

Tlj" ■ 1 Introduction 

o 

:^ 

^ • Heavy quarks are those with masses ruh S> Aqcd such that as{rn\) oc 

■ ln~"'"( ) ^ 1 and hence, according to this definition, the charm, bottom and 

1^ ' top quarks (c, 6, t) are heavy whereas the up, down and strange quarks (m, d, s) are 

' hght. Apart from the fact that heavy quark production processes are fundamental 

elementary particle processes which take place with substantial rates at high cn- 
k>( \ ergy colliders, and which therefore require a good phenomenological understanding, 

■ heavy quark production is theoretically interesting for at least two more reasons. 
I (i) Firstly, the heavy quark mass ruh acts as a physical long-distance cut-off such 

that the heavy quark (or more precisely logarithms of the heavy quark mass) can be 
treated within perturbative QCD - either within fixed order perturbation theory or 
in resummed approaches where a heavy quark parton distribution function (PDF) 
is introduced using perturbatively calculated boundary conditions for the Altarelli- 
Parisi evolution equations, (ii) Secondly, in addition to ruh another large scale, e.g. 
the transverse momentum pT of the heavy quark, is usually involved in the hard 
scattering process such that one has to deal with multi-scale processes which com- 
plicate the perturbative analysis. Obviously, terms of the order 0{m\/p'^) should 
not be neglected in resummed approaches if pt is not much larger than the heavy 
quark mass. This remains also true if the perturbative boundary conditions in (i) 
would be abandoned in favour of experimentally determined input distributions for 
the heavy quark PDF. 

The fixed order treatment is also called fixed flavour number scheme (FFNS) 
since the number of flavours in the initial state is fixed to Uf = 3(4) for charm 
(bottom) production. On the other hand, in the parton model the number of 
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active flavours is increased by one unit, Uf + if the factorization scale 

crosses certain transition scales (which are usually taken to be of the order of the 
heavy quark mass)". Accordingly, these schemes are called variable flavour number 
schemes (VFNS). The conventional massless parton model is usually referred to as 
massless or zero mass VFNS (ZM-VFNS) whereas massive parton model approaches 
are named massive or general mass VFNS (GM-VFNS). For details see, e.g., Refs. 

I2EI- 

In this contribution I will consider open charm hadroproduction, pp — > D*X, 
and describe a new calculation in a massive VFNS |3] guided by the factorization 
theorem of John Collins including heavy quark masses 0. This work is along the 
lines of previous studies of the direct |6j and the single-resolved contributions |7] to 
the process 77 — > D*X. Note also that the latter constitutes the direct contribution 
to 7p D*X considered in |8I3| . The results of this new calculation can then be 
used also to obtain the double resolved contribution in the 77 and the resolved 
contribution in the 77? process in order to perform a complete NLO analysis of 
these processes in this massive VFNS. 

Recently, the CDF collaboration has published first cross section data for single 
inclusive D meson production in pp collisions [H] obtained in Run II at the Tevatron 
at center-of-mass energies of ^fS = 1.96 TeV. These data have been compared with 
calculations |HJI11| in two distinct theoretical approaches, the FONLL approach ^21 
and the conventional massless parton model which both are compatible with 
the data within theoretical and experimental errors.'' The work presented here 0] 
extends the latter successful calculation in the massless approach by terms of the 
order 0{m1/p^). 



2 Basic Fi"amework 

Our theoretical basis for the calculation of differential cross sections for open charm 
hadroproduction, pp D*X, is the familiar factorization formula [see Fig. 
however, with heavy quark mass terms included in the hard scattering cross sections 



"For a detailed discussion see the appendix in [J and references given there. 

''However, for the central values the ratio of data/theory is about 2 at low pj- and 1.5 at high pj-. 
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Here ff{xi,iif) and /j'(a;2, /i/), i, j = u,d,s,c,g, are universal PDFs of the proton 
and the anti-proton which are non-perturbative quantities in the case of gluons, 
hj = ffj Slid hght quarks, i,j = q = u,d,s. On the other hand, the charm 
quark PDF is usually obtained from perturbatively calculable boundary condi- 
tions. Furthermore, the fragmentation of the final state parton k = g,q,c, . . . 
into the D* meson is described by universal non-perturbative fragmentation func- 
tions (FFs) {z,fj.'j:). Finally, the hard scattering cross sections d&{ij 
kX){iJ,f , fi'j:, as{fi), ^) are perturbatively calculable and depend on the factoriza- 
tion scales (If and /z^ and the renormalization scale ^. As has been shown in Ref. 
[5], heavy quark mass effects can be consistently included in the hard part where 
it is of course mandatory to use the same factorization scheme for PDFs, FFs and 
the hard part. 

3 Hard scattering coefHcients with masses 

For the calculation of massive hard scattering coefficients we have adopted the 
following procedure |6I7I3| : (i) We have calculated the m limit of the massive 
3-FFNS coefficients |14I15| only keeping m as regulator in logarithms In ^ . The 
partonic subprocesses occurring in the 3-FFNS are gg — > cc and qq cc in LO 
and gg — s- ccg, qq — > ccg, and gq ccq (gq — > ccq) in NLO. Here special care was 
required in order to recover certain distributions ((5(1 — w), ( '"[1_~^"'* ) + ) 

occurring in the massless MS calculation, (ii) Subsequently, we have compared 
the massless limit with the corresponding short distance coefficients in the genuine 
massless MS calculation in order to identify the subtraction terms by the unique 
prescription 

dcrsuB = ( fim dcr(m)) - d<T(MS) . (2) 

rn — *0 

The massive hard scattering cross sections have then been obtained by removing 
the subtraction term from the massive coefficient in the 3-FFNS: 

da{m) = da{m) — dcsuB ■ (3) 

It should be noted that by this procedure no finite mass terms are removed from 
the hard part apart from the coUinear logarithms In ^ . (iii) Contributions with 
charm quarks in the initial state have been included in the massless approach |15j . 
Note that this can be done without loss of precision ^^l- Moreover, this rule is of 
great practical importance since the existing massless results for the hard scattering 
coefficients |17| can simply be used whereas their massive analogues are unknown 
and would require a dedicated calculation of these processes. Note that massive 
heavy quark initiated coefficients have been obtained for the case of deep inelastic 
scattering jl8ll9| and the results are already quite involved. Finally, this rule 
(scheme) has also been employed in a recent analysis of CTEQ PDFs with heavy 
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quark mass effects |5U] such that it is necessary to do the same if these PDFs are 
to be used in a consistent analysis in the future. 

4 Numerical Results 

In this section numerical results for the gg ^ c + X channel (i.e. for pp — s- {gg 
c) D* + X) are presented in form of the ratio 

^ _ I-i'iy dggg/dydpT 

ill dy dal^Q/dydpTim = 0) ' 

in dependence of the pt of the meson. In order to investigate the size of the 
mass effects the cross section da/dpTdy in the numerator has been calculated in 
two ways: firstly, in the massive 3-FFNS jTJ using the FORTRAN code provided 
by I. Bojak and secondly, using the massless limit which we derived analytically 
as described in the preceding section. In both cases the LO cross section in the 
denominator has been calculated with massless charm quarks. 




Figure 2. Contribution of the gg — > cX channel to pp D* + X normahzed to the LO cross 
section with m = 0. Shown are results for the massless limit (solid line), the massive calculation 
(dashed line), the subtracted massless limit (dotted line), and the subtracted massive calculation 
(dashed-dotted line). The Factorization and renormalization scales are (a) fif = fi'j: = /ir = 
(but fixed at /i = 2.1m in PDFs, FFs, and as) and (b) /if = fi'r = fir = 2my. 



The results in Fig.|3have been obtained using the CTEQ6M PDFs and the 
FFs for c D*^ from (22j (NLO OPAL version) along with the following input 
parameters, m = 1.5 GeV, Ep = Vs/2 = 980 GeV, A(„^=4) = 328 MeV (i.e., 
cts{mz) = 0.118). In Fig. |21(a), the renormalization and the factorization scales 
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are fir = fJ'f = fJ''f = itt-, but the scale in as, in the PDFs and the FFs is chosen as 
/i = 2.1m, to stay above the starting scale in the PDFs and FFs. The solid lines are 
the massless limit and the dashed curves are the result of the massive calculation. 
As one can see, the massive cross section approaches the massless one very slowly 
at large pt- At px = 20 GeV the difference between the massive and the massless 
result is still of the order of 6% and the massive cross section is always larger than 
the massless cross section ratio in the pr range between 7 and 100 GeV. From this 
comparison we conclude that in the gg channel the terms proportional to m? /p^ are 
particularly large and lead to an increase of the massive cross section as compared 
to the massless approximation. This agrees with findings in jl2| , where the massive 
and the massless theory were compared as a function of the mass m for fixed pt- For 
our application we are interested in the massive cross section and its massless limit, 
where the subtraction terms dcr^^^ are removed. This leads to the dashed-dotted 
(massive) and the dotted (massless) curves in Fig. 12 (a) . Not shown is the genuine 
MS result, generated with the routine of |17I23I24| . which is in exact agreement 
with the subtracted massless limit (dotted curve) demonstrating the correctness of 
the subtraction terms for /i^ = A'/ — M/ — Fig- 121(b) shows results for the choice 
fir — ^J■f — ^J''f — 'iviT = 1^ w? + pip which can be used for small and large pt- 
The cross sections for these scales are shown again in four curves, massless theory, 
massive 3-FFNS, subtracted massless theory and subtracted massive theory where 
the labeling of the curves is the same as in Fig. |21(a). These curves are our final 
results for the gg channel. As can be seen, for this different scale the cross section 
ratios are much flatter as a function of pt than in Fig. [3 (a). The K factor is 
somewhat smaller now, but it is still large, showing that for the gg channel the 
perturbative expansion is not converging very well. As in Fig. |2] (a) we observe 
that the massive cross section converges from above to the corresponding massless 
cross sections with increasing pT ■ The ratios for the subtracted cross sections are 
now negative for px < 30 GeV and rise up to « 0.6 at pT = 100 GeV. Also for these 
scales we find exact agreement between the genuine MS result and the subtracted 
massless limit (dotted curve). Thus we can be sure that we have subtracted from 
the massive cross section the correct terms in order to get the massive VFNS cross 
section in the MS subtraction scheme, which approaches the massless limit in this 
scheme at large pr- 

In contrast to the gg channel, the mass effects in the qq and gq channels are 
quite small reducing the size of the mass effects in the sum of all three channels. 
Furthermore, contributions with charm quarks in the initial state and contributions 
from the gluon (and light quarks) fragmenting into the D* meson have to be in- 
cluded. These contributions are dominanf^ and calculated with m = as explained 
in Sec. 13 This further dilutes the mass effects occurring in the gg channel such 
that quite small and positive mass effects will survive in the total result 4 . 



■^This can be understood already with help of Fig.|2](b) where large subtractions are visible which 
are resummed in the evolved charm quark PDFs (and FFs). 
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